Globular Structures of a Helix-Coil Copolymer: 
Self-Consistent Treatment. 
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A self-consistent field theory was developed in the grand-canonical ensemble formulation to study 
transitions in a helix-coil multiblock globule. Helical and coil parts are treated as stiff rods and 
self-avoiding walks of variable lengths correspondingly. The resulting field-theory takes, in addition 
to the conventional Zimm-Bragg parameters, also three-dimensional interaction terms into account. 
The appropriate differential equations which determine the self-consistent fields were solved numer- 
ically with finite element method. Three different phase states are found: open chain, amorphous 
globule and nematic liquid-crystalline (LC) globule. The LC-globule formation is driven by the 
interplay between the hydrophobic helical segments attraction and the anisotropic globule surface 
energy of an entropic nature. The full phase diagram of the helix-coil copolymer was calculated and 
thoroughly discussed. The suggested theory shows a clear interplay between secondary and tertiary 
£N| structures in globular homopolypeptides. 

PACS numbers: 61.30.Vx Polymer liquid crystals; 87.14.Ee Proteins; 87.15.-v Biopolymers: structure and 
M-H physical properties 
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I. INTRODUCTION 



Significant work has been done to investigate the helix-coil transition theoretically-^ and computationally^^. One 
of the first and well-known approaches is the Zimm-Bragg theory^ which is designed for single-stranded polypeptide 
chains. It considers a one-dimensional Ising type model in which each segment can be in one of the two states: helical 
state (stiff) or coil state (flexible). The helix is stabilized by hydrogen bonds which generates an energy gain of — e 
for each segment in the helical state. This energy gain is partly balanced by an entropy loss —AS. The free energy 
difference between the helical state and the coil state is therefore given by A/ = — e + TAS for each segment. 

In an ev-helix hydrogen bonds can be only formed between the i th and the (i + 3) th peptide group. The formation 
of a hydrogen bond between the first and the third peptide group fixes the conformation of three groups. The next 
bond between the second and the fourth group furnishes the same energy gain but fixes only one new group and thus 
\Q leads to a much smaller entropy loss. The formation of an a-helix ( as a simple element of the secondary structure) 
is therefore a cooperative process and the formation of a junction between helix and coil is energetically unfavorable. 
This can be modelled by an energy penalty /ij for each junction between two successive groups of helical and flexible 
segments. Similar arguments hold for all kind of helices, such that the helix formation is always a cooperative process. 
It is convenient to define the following fugacities 



s = e-M, a = e- 2 ^\ (1.1) 



where s gives the statistical weight of a helical segment compared to a coil segment. The cooperativity parameter 
a gives the statistical weight of a junction point, a = 1 corresponds to /xj = and therefore to a non-cooperative 
system, a — > corresponds to [ij — > oo and therefore to a totally cooperative system, i.e. either the entire chain 
forms one big helix or no helix is formed at all. In most helix forming biopolymers a is roughly 1CP 3 — 10~ 4 (see e.g. 

Using a transfer matrix method^ the one dimensional model can be solved exactly. The s - dependance of the 
fraction of stiff helical segment Or = N-r/N has a typical sigmoid - type form . With increasing cooperativity (i.e. as 
the parameter a is reduced) the transition becomes sharper and sharper. This rather sharp crossover transition (due 
to the cooperativity effect) is also observed experimentally, for instance in polybenzylglutamatei 

Several extensions have been made to this one-dimensional Ising type model. It has been shown^ that the transition 
becomes less cooperative, if the hydrogen-bonding ability of the solvent is taken into account. The helix-coil transition 
in grafted chains was studied^ as well as the effect of an external applied force on the transition^. 

Of special interest is the application of one-dimensional models to proteins, see for instance^. A study of the helix- 
coil transition including long-range electrostatic interactions^ can, to some extent, explain the amount and location 
of helical segments in globular proteins. However, to understand how cv-helices (or generally secondary structure 
elements) are formed in the folding process of a protein and how this influences the compaction and formation of 
tertiary structure (and vice versa), it is necessary to combine the one-dimensional physics of the helix-coil transition 
with the interactions of segments which approach each other in the three-dimensional space. This enables a description 
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of the interplay between the secondary structure and the mesoscopic structure formation of the entire chain (known 
as tertiary structure) . 

In proteins the stiff helical parts are often hydrophobic since the hydrogen bonds stabilizing the helix composition 
(as opposed to coil parts) are mainly saturated (due to bond formation between the i th and the (i + 3) th peptide 
groups). This hydrophobicity causes an additional attractive interaction between stiff parts and drives the protein 
into a compact globular phase. Statistical analysis of the data from 41 globular proteins in native and partially 
folded conformational states^ showed a strong correlation between the amount of secondary structure elements and 
compactness of the proteins. This indicates that the formation of secondary structure (for instance a-helices) and 
the hydrophobic collapse into a compact globule occur simultaneously. This problem has been partially discussed 
within computer simulations of globular protein a 12 i 13 ' 14 i 15 . Among other things, it has been show n 13 ' 14 i 15 (as opposed 
to the earlier findings by Dill et alJ^) that the compactness itself, driven by the hydrophobicity, is insufficient to 
generate any appreciable secondary structure. It was necessary to introduce a local conformational propensity toward 
a - helix formation. The interplay of compaction and secondary structure leads to the formation of the specific 
three-dimensional tertiary structure. Computational and experimental studies of this mechanism can for instance be 
found in ref i 16 ' 17 . 

To study the interplay of helix-coil (or stiff-flexible) transition and collapse transition of the polymer into a compact 
globule, we have developed an approach which combines variable composition with three-dimensional excluded volume 
interactions using self-consistent field theory. The paper is organized as follows. In Section II we have covered the 
general self - consistent field theory of a helix - coil copolymeric globule. The final equations are written down in a form 
which is convenient for the numerical solution. Section III is devoted to the analysis of this solution. The formation 
of different globule structures (e.g. amorphous and liquid - crystalline (LC) globules) are studied in details. Among 
other things we argue that the LC - globule formation is mainly driven by the globule surface tension anisotropy. 
Finally, in Section IV we sketch the main results and compare them with the appropriate findings of some other 
authors. 



II. THEORY 

In this section the derivation of the field theory is outlined. For a more detailed description see 1 ^. Regarding all 
three dimensional interactions and entropic contributions the helices are modelled as stiff rods. Hence for one specific 
microscopic composition of helical and flexible parts the system looks like a rod-coil multiblock copolymer as shown in 
Fig-©- The multiblock copolymer may be composed of K rod-coil blocks. The conformation of the n th rod-coil block 
is given by the vector- function r„(s) describing the contour of the coil, by the vector R„ which gives the position of 
the junction point between rod and coil and the unit vector u„ describes the orientation of the rod, see Fig.©. The 



r„(f„NJ 




FIG. 1: Rod-coil multiblock copolymer parameterization. Two successive rod-coil blocks with ordinal numbers n — 1 and n are 
shown. The flexible chain coordinates use the vector set r(s). The orientation of the n th rod is denoted as u n and the junction 
point between rod and coil is given by R„. 

length of the n th rod-coil block in units of the segment length b is given by N n . The fraction of the flexible segments 
(coil) is given by f n and the fraction of stiff segments (rod) thus by 1 — /„ . The microscopic flexible segment density 
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pc( r ) an d stiff segment density p R (r) can now be defined as follows. 

K f«N n 



PC 



(r) = J2 I ds5(r-r n (s)), 

n=0 { 



(2.1) 



p R (r) = ^2 / dsS(r - R n - u n sb). 

M — n » 



(2.2) 



In addition an orientation density S lJ (r) is introduced which is sensitive to the collective orientation of the system 

(l-U)N n 

(2.3) 



S« (r) = / «J(r - R„ - u„s6) 

n 



n=0 



The interaction Hamiltonian of the model can thus be written as 



-S ij 



f3H int [p c , PR ] = X J d 3 rp R (r)p R (r) + ^ J d 3 r[ Pc {r)+p R {r)f 

+ i / dSr [pc(r)+PRir)f + 9 f <^Tr [s«(r)S«(r) 



(2.4) 



where v and w control the strength of the global two- and three-body interactions. The parameter \ controls a 
selective two-body interaction between the stiff segments which is caused by the hydrophobic nature of the helical 
parts of the chain (cf. Introduction). The last term is a alignment interaction between the rods of the standard 
Maier-Saupe fornJ£. The representation of this type has a wide use in the polymeric liquid-crystal contexfj^S. 

The canonical partition function of the entire system can be written as the functional integral over the collective 
coil, rod and orientation densities, pc{ r ), Pr{. v ) and S tJ (r) respectively, i.e. 

Z({N n },K) = fv Pc {v) fvp R {v) fvS ij (r) f ]]_ Dr n (s) d 3 R n d 2 u n 
J J J J n _ 1 

x S(p c (r) - pc(v))6(pr(v) - fe(r))5(5* J '(r) - &*(t)) 



!f n N n , 
"i / ds {it) 


x exp I -x / d 3 r p R (r)p R (r) - ~ / d 3 r [p c (r) + p R {r)} 2 



w 
3! 



d 3 r [p c (r) + p R (r)} 3 - g / d 3 r Tr [S«(r)S«(r)] 



(2.5) 



By making use of the integral representation for the S - function in Eq. (|2.5[) (which results in an appearance of external 
fields hc{r), h R {r) and h l g (r)) the partition function of the entire system reads 

Z({N n },K) = J V Pc (r) Jv Pr {y) J DS^(r) J Vh c (r) J Vh R (r) J Dh^r) 

x exp|-x J d 3 rp R (r)p R {r)-^ J d 3 r [p c {r) + p R {r)} 2 

d 3 r[p c (r)+p R (r)} 3 -g J d 3 r Tr [S ij (r)S ij (r)] 

+ i J d 3 rp c (r)h c (r)+i J d 3 r p R (r)h R (r) + i J d 3 r S ij (r)/$ (r)| 

x Z^({N n },K-[h c },[h R },\h^}) , (2.6) 
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where the partition function of the non - interacting system in the external fields he, /ir, h l £ . is given by 



Z^({N n },K;[h c ],[h H ], 4 ) 



K 



x exp 



Y[Vr n {s)d 3 R n d 2 u n S(\u 

71=1 

3 

262 



fnN n 



ds 



ds 



l)S(r n (f n N n )-H n ) 
dsh c {y n {s)) 



ds /ir(R„ + u„s) 



(l-f„)N n 



ds h l g (R„ + u. n 



« 



(2.7) 



The composition of the system is assumed to be equilibrated with respect to the total number of stiff (Nn) and 
flexible (Nc) segments as well as to the number of junction points (-/Vj) between stiff and flexible parts . Therefore 
the description can be reduced, i.e. {N n }, K — > TVr, Nc, Nj, and it is more convenient to switch to a grand canonical 
partition function Z of the interacting polymer written in terms of the grand canonical partition function Z^ of the 
non-interacting polymer in the external fields he, /ir, h l g , i.e. 



Z(n,e,tr) = J V Pc (r) J Vp R (r) J DS ij (r) J Vh c (r) J Vh R (r) [ Dh%(r) 



x exp \-pH mt \p c ,p K ]+i / d\ Pc {v)h c {v)+i / d 3 rp K (r)h R (r)+i / d 3 r^(r)^'(r) 



x ZW((,,e,a;[h C ],[h K ], ft« ). 



(2.8) 



In Eq. (l2.8j) p denotes the chemical potential conjugated to the whole number of segments and — e is the energy gain 
of a helical segments compared to coil segment. The meaning of the cooperativity parameter a was discussed in the 
introduction. 

The grand canonical partition function Z^> of the non-interacting system can be derived by using the polymeric 

) = / dldl / a(°)(l > l';/i > e > ff;[/»c],['iR]' 



correlation function 3® , i.e. Z^ (fi, e, a; [he] , [hn] , 
meric correlation function E^ -* (1, 1'; p,, e, a; [he] , [fin] , 



h« 

'4 



h '■> 



The poly- 



) gives the conditional unnormalized probability of finding 



the first segment of the copolymer at the coordinate 1 provided that the last segment is at 1'. The symbol 1 stands 
either for ri or for (ri,Ui) depending on whether the first segment is a flexible or a stiff one. The same holds for 
the coordinate 1' of the last segment. The polymeric correlation function is therefore the partition function of the 
multiblock helix-coil copolymer with two ends fixed at 1 and 1'. 

By using the equations of motion for the rod and coil Green functions - G co ii and G roc i - the polymeric correlation 
function S^ ^ can be constructed. The inverse Greens function operators are given by^- 



G-i, = 5(r-v')((3p- 



b 2 . 



G7ad = 6(r-r')(p(p~e)~ 



ih c {r) 



& 2 (u-V R ) s 

(3(p - e) 



+ ih K (r) + ih s (r) : P ), 



(2.9) 
(2.10) 



where the tensor P is given by P^ 



S tJ /3. In Eqs. (|2.9p - (|2. 10[) we have set the elementary unit lengths in 



coil and rod parts equal to each other (namely to b) for simplicity. This affects some quantitative values but definitely 
does not alter qualitative predictions of this paper. 

Knowing the inverse Green function operators the grand canonical polymeric correlation function can now be 
represented as the following Gaussian 2-dimensional path integral 



S(°)(l,T;/i,e,a; [h c 



h l i 



) = i / P,-ZVr(1 )v il') 



x exp ■ 



d3d4 



^(3) 
U( 3 ) 



^rod 

__l/2 (~t-\ 
~° ^coil 




(2.11) 
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where 



6 = / VipVipexp < -- 



1 



d3 d4 



^(3) 
U( 3 ) 



^rod 
^.1/2 




(2.12) 



The choice of -0(1) vO-') under the path integral yields one rod-coil unit as a basic building block, see Fig.® and 
Eq. (|2.14p . This choice assigns the coordinates (1) and (1') to (r, u) and (r') respectively. Similarly, and ^(4) 
are shorthand notations for ^(r 3 , u 3 ) and ip(r 4 , u 4 ). 
The inversion of the 2 x 2-matrix in Eq. (|2.1ip reads 



f}- 1 



rV2 



a ^coil 



^rod * ^coil 



a 1 / 2 G-\ 



(2.13) 



The calculation of the path integral in Eq. (|2.11[) yields the following result 



sW(l,l';M,c,ff; [he], [ha], 



^rod * ^coil - ° 
= cr 1 / 2 G r od * G co ii * 



1 + crGrod * G C oil 



f G ro d * G co il * G ro d * G C oil + 



(2.14) 



The asterisk in Eq. (|2.14|) is a shorthand notation for a convolution of Green function operators. The geometric 
progression with a convolution as a binary relation has a clear pictorial representation, see Fig. ((21). The first term of 



(l, l';[h c ],[h R ],[h]) = 



+ 



+ 



FIG. 2: Pictorial representation of the geometric progression in Eq. (|2.14p . One rod-coil unit constitutes the basic building 
block. 



this series has the analytical expression cr 1 / 2 G ro d * G co u, whereas the ratio is equal to crG ro d * G co n, see also Fig. ([3]). 
The first term represents one rod-coil unit with one junction between rod and coil, hence the factor a 1 / 2 . If an 
additional building block is added, two more junctions are created, hence the factor a in the ratio. This series gives 
a correct representation of the grand canonical polymeric correlation function. 




A 



A 



A 

G 



= a G rod * G 



coil 



coil 



FIG. 3: Pictorial representation of the series ratio in Eq. (|2.14[l : each bar (rod), zigzag line (coil) and fat dot correspond to 
Grod,G co ii , and a 1 ^ 2 respectively. 
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The denominator in Eq. (|2.1ip can be avoided by introducing de Gennes' n — > trick 21 . Consider the two n- 
component vector fields {ipa, fa}, where a — 1, 2, ... n. Then Eq. (|2.11[) can be formally rewritten as 



n „ 

H(°)(l,l'; M ,e,a; [he] ,[hn] , [h%] ) = lim J] / Thl> a V<p a ^(l)^i(l') 



~ a °coil 



Va(4) 



(2.15) 



Integrations over the external fields he, hn,, h l g, over the densities pc(r), /OR(r), S 4J '(r) and over all endpoints {1, 1'} 
yield the final field theoretic representation of the grand canonical partition function 

Z(n,e,a) = Urn f[ J V^ a Vip a J d 3 r d 2 u ^(r, u) ^ dV'^^r') 



exp /d 3 rd 2 «i/> Q (r,u) /3{/J.-e) 



b 2 (u • V r )' 



- e) 



V> a (r,u) 



/ d 3 ^ a (r)k--vi„(r) 



~ 7 /WE /d 2 «^(r,u) 

/Ti n II 

d 3 r £ /d 2 W ^(r,u)+^^(r) 

r n 3 

- S/ d ' p E/d 2 U V^(r,u)+X:^(r) 

La— 1 a' — l 

n « 

+ ^ 1/2 E / d 3 rd 2 M Va(r,u)^ a (r) 

Q = l ^ 

n n p ~\ 

d 3 rd 2 W dVF 2 (u.u')C(r,u)^(r,u') 

a=ln'=r J 
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where the second Legendre polynomial is given by 



P 2 (0) = -(3cos 2 0-l) 



(2.16) 



(2.17) 



To evaluate the partition function Z(/j,, e, a) the self-consistent field approximation is used22. In the self-consistent 
field approximation fluctuations are neglected and the functional integral over the fields in Eq. (|2.16p is integrated by 
steepest descent. The saddle point solutions for </5«(r) and i/j a {r, u) are chosen such that the effective grand potential 
keeps the full symmetry of the Hamiltonian in replica space, i.e. it is invariant under rotations in replica space. This 
is the case for (see ref^i) 



V>a(r,u) = n a ip(r,u) 
tp a (r) = n a <p(r), 



(2.18) 



where n is a unit vector such that n a — 1- 

In order to make the problem more tractable we expand tjj(r, u) in spherical harmonics, ip( r , u ) = 
m ipi m (r)Y[ m (u.), (see e.g. ref. 23 ). Since the solution for i/?(r,u) must respect uniaxial and cylindrical symme- 
try, the expansion reduces to Legendre polynomials, i.e. m = 0. This expansion is truncated to the lowest nontrivial 
order— 



1 



4tt 



1/2 



^(r,u)« — M*)+ hr V- 2 (r)P 2 (u-n z ) 



4tt 



1/2 



(2.19) 
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where n is the main direction along which the rods in the core of the globule are aligned, if the system forms an 
anisotropic globule with aligned rods. If they are not aligned and the system forms an amorphous globule, only the first 
term in the expansion in Eq. (|2.19p differs from zero. On the other hand, the second term in Eq. (|2.19p is responsible 
for the nematic LC-order. The main direction of alignment n can be chosen arbitrarily without loss of generality, since 
a change in alignment direction only corresponds to a rotation of the complete globule in the laboratory coordinate 
frame. Thus we choose n directed along the z-axis of the (x,y,z) laboratory frame, i.e. n = n z . This expansion 
allows us to perform the u-integrations. 

The resulting effective saddle point grand potential il(/i, e, cr) can now be calculated. It is given in terms of the 
saddle point fields tpo(r), ^2 (r) and (p (r) by 



/?(/! - e) 



d 3 r [V> 2 (r)+V 2 2 (r)] 



210/3(jU-e) 



d 3 r{35^ (r)V 2 Vo(r) 



2 di -di- di 



Mr) [25 d 2 x + 25 8 2 y + 55 3 2 ] V 2 (r)} 



+ 



d 3 r^(r) pn- -V 2 ip(r) 
L J 



w 
48 



d 3 r [^ 2 (r)+^(r)] z 
d 3 r [^ 2 (r)+^ 2 (r)+^ 2 (r)] 2 
d 3 r [^(r)+^(r)+^ 2 (r)] 3 
d 3 rip(r) ip (r) 

+ 245 / d3r {^W [7^o(r) + V5^ 2 (r)]} 2 . 

The coil and rod densities as well as the orientation density are given by the following relations 

1 



(2.20) 



Pc(r) 
S zz (r) 



\ J ^ 2 (r,u),i[f 3 (r 

i J d 2 iiP 2 (cos0)V 2 (r, u) 
2 r 

—7=ip2(r) Vo(r) +V5ip 2 ( r ) 
V5 L 



^ 2 2 (r)] 



(2.21) 



Functional minimization of Eq. (|2.20p with respect to ipo(r), ^2(r) and ip(r) yields a set of three coupled partial 
differential equations, which are solved numerically with the finite element tool kit Gascoignc 25 . The results are 
extensively discussed in the next section. 



III. RESULTS 



In this section the numerical results for the full set of equations describing a rod-coil multiblock copolymer with 
a variable composition of stiff and flexible segments are presented. The segment length is set to b = 1. In addition, 
all energies such as e, /1 and also the saddle point free energy F are given in units of k^T . In this section this will 
not be indicated in order to avoid complicated notation. Eq. (|2.20[) describes the copolymer in a grand-canonical 
representation. In the grand-canonical ensemble the total number of segments of the polymer N is not fixed but its 
mean value is determined by equilibrium conditions. 

A real polymer has a fixed length. In order to ensure this fixed length TV, the chemical potential \i is - for each set 
of physical parameters (v, w, x, g, e, cr) - tuned such that the equilibrium value of TV is equal to the desired one. The 
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total number of segments N — N co u + N T0 a is calculated by numerical integration over the rod and coil densities given 
by Eq. (|2.21[ ) . For a given set of parameters, N(/i) can be computed 2 ^ and a typical example of this curve is shown 
in Fig. ([4]). For \i — > the total number of segments N diverges. This corresponds to the N ~ behavior which is 
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FIG. 4: N as a function of p, for w = 1, v = -0.2 and a = 10" 4 . The dotted curve corresponds to e = 0.08 and \ — 0- The 
continuous curve corresponds to e = 0.1 and \ — 0.0138. 



well-known for a O-solvent chain 2 ^. The divergence of N at a specific value /z on the right hand side of the minimum 
corresponds to a fully collapsed infinite globule and has been discussed first by Kholodenko and Freed2i. These two 
branches of iV(//) meet each other in the minimum of N(/j,) which can be associated with the coil-globule transition 
point, as will be shown in the next section. Since N is always fixed by tuning /i, it is possible to distinguish from a 
plot like the one shown in Fig.(0| whether the system is left of the coil-globule transition point (i.e. in the open chain 
regime) or right of the transition point (i.e. in the globular regime). However, for fixed N, it is necessary to choose 
one of the two branches. Since this work focuses on the study of globular structures, the numerical calculations are 
always restricted to the right branch including the minimum. The self-consistent field theory is only expected to give 
good results for this branch since fluctuations are neglected. 

The three-body interaction parameter w is chosen to be w — 1 throughout the entire paper and the two-body 
interaction constant v is always negative to ensure that the system stays in the globular regime (up to the transition 
point). 

As a first step the pure homopolymer globule case^i has been studied in order to test the numerical routines. For 
this end one should simply set = "02 = in Eq. (|2.20p . Here we are not going to discuss these results in detail . 
It is pertinent only to note that the resulting numerical solution fully supports the following well known theoretical 
results. The critical value of the two-body interaction constant, v, scales as |u c | ~ A -1 / 2 , whereas the maximal 
globule density in the critical point behaves as p cr it ~ N~ 1 ' 2 . 

A short remark on the terminology that will be used below is necessary at this point. The terms phase and transition 
will be used frequently although the system is a polymer of finite length. All transitions are therefore crossovers of 
finite width with continuous order parameters. It is nevertheless common now as applied to "soft matter" to use the 
the term "phase" to distinguish the different structural states of a polymer and to refer to the crossover between these 
states as a "transition". 



A. Coil-globule transition 



The rod-coil multiblock copolymer shows a coil-globule transition similar to the one of a homopolymer. To demon- 
strate this, the interactions which are specific for the stiff segments are set to zero, that is x = 9 = 0. In addition, the 
energy gain per stiff segment is set to zero (e = 0) and it is assumed that there is no cooperativity in the formation of 
stiff segments (a = 1). The two-body interaction constant v is varied. The transition point between coil and globule 
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is defined to be the minimum of the iV^-curve as it was explained above. The length of the polymer is fixed at 
N = 550. Fig. ([5|) shows the profile of the total density p(r) of the copolymer as a function of radial distance from the 
center. The total density p(r) at each point is given by 



p(r) = pc(r) + p R (r) = ^ 2 (r) + ~^(r) + ~V>f(r). 



(3.1) 



As can be seen from Fig.©, the density profile becomes broader with decreasing \v\. At v = —0.5 the copolymer is 
1 




V — 


-0.5 




-0.3 




-0.077 



20 



25 



FIG. 5: This plot shows the radial density profile of the entire copolymer in radial direction for different values of v. The 
dashed curve for v = —0.077 corresponds to the coil-globule transition point. 



deep in the globular state with a big plateau of almost constant density and a rather small surface layer of decreasing 
density. At the transition point v — —0.077 the plateau basically vanished and the surface layer becomes very broad. 
To further illustrate the structural change, Fig.© shows a colour-coded plot of the local density in g-z space, where 
g denotes the radial direction and z the axial direction in cylindrical coordinates. The center of the globule is located 
at the bottom left corner of each picture. The pictures on the left show the local density of flexible segments and the 
pictures on the right correspond to the local density of stiff segments. Red indicates high density and dark blue zero 
density. In the upper two plots the copolymer is at the transition point (v = —0.077). In the lower two plots it is 
deep in the globular state (v = —0.5). 

A much clearer indication that v = —0.077 corresponds indeed to a transition point can be seen from Fig. ([7]), where 
the globule radius is plotted versus v. The globule radius is defined as the point R in radial direction at which the 
density p(R) has decreased to p(R) = 10 _3 po, where po is the maximum density at the center of the globule. The 
radius R shows a rapid increase when v = —0.077 is approached. Note, that the copolymer is finite (here: A=550) 
and therefore all transitions are crossovers as discussed above. 



B. Helix-coil transition 



In this section the fraction of stiff segments is investigated as a function of the energy gain per stiff segment e. 
This is similar to the helix-coil transition described by the Zimm-Bragg model 4 . A major difference is that the 
model used here is a three-dimensional model of the polymer including interactions, whereas the Zimm-Bragg model 
and its extensions are one-dimensional models with no three-dimensional interactions and no explicit entropy term 
. On the other hand, the Zimm-Bragg model can be solved exactly, whilst the self-consistent field treatment of the 
three-dimensional model is a mean-field approach which neglects fluctuations. 

Two different regimes will be discussed in the following: a low cooperativity regime with a in the range 0.05 — 1 
and a high cooperativity regime with a in the range 7 ■ 10~ 3 — 10~ 4 . Remember, that a — 1 means no cooperativity 
and a — means total cooperativity. Throughout this section \ and g are set to zero. There are therefore no specific 



10 




FIG. 6: The density of flexible and stiff (helical) segments are shown on the left and on the right panels respectively, v — —0.077 
corresponds to the upper plots and v — —0.5 to the lower plots. 
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FIG. 7: The globule radius R is plotted as a function of v. 



interactions between the stiff segments. The only interactions are attractive two-body interactions and repulsive 
three-body interactions between all segments. Fig.© shows how the fraction of stiff segments Or depends on e for 
different values of a in the range 0.05 — 1. Even for small cooperativity, the slope clearly depends on a and gets larger 
with increasing cooperativity (decreasing a). 
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The high cooperativity regime is shown in Fig.©. The slope of the 0R,(e)-curves increases with decreasing a as 
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FIG. 9: The fraction of stiff segments 0r is plotted as a function of e for high cooperativity. v = —0.2 and N = 10 5 . 

expected. Especially for a — 1CP 3 (green curve) it can be seen that the curve is asymptotically approaching 0r = 1. 
The curves for a = 10~ 4 , 5 • 10 -4 and 10 -3 end at a certain value of e. These values correspond to the coil-globule 
transition point. This can be explained as follows. If there is no additional selective interaction energy which favors 
a compactification of the stiff segments (i.e. x — .9 = 0), the stiffening of parts of the chain due to an increase of 
Or with increasing e pushes the chain segments further apart from each other and therefore leads to a more open 
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structure. From a certain value of e on, the system is thus pushed into the open chain regime. For higher cooperativity 
this effect is stronger, since the system forms less junctions points and hence on average longer rods. For smaller a 
the coil-globule transition point is therefore reached at smaller values of Or and e. 

To illustrate this stiffening, in Fig. (fTU|) the average rod length Lr is plotted as a function of e for a — 1CP 4 . In 
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FIG. 10: The average rod length Lr is plotted as a function of e for a = 10" 4 , v = -0.2 and TV = 10 5 . 



units of b the average rod length Lr is given by the total number of stiff segments divided by half the total number 
of junction points between stiff rod and flexible chain. In the self-consistent field theory approach Lr is given by 

Lr = l£dggj^dz [y>g(g, z) + tgjg, zj\ 
b f °° dggf^dz <p(g, z) ip (g, z) 

Fig, (fit))) shows that for e = the average rod length is roughly equal to 1. For e > the cooperativity effect sets 
in and the average rod length increases up to Lr » 2.3. At this value the stiffening of parts of the chain is strong 
enough to drive the polymer in the open chain regime. 

It is of interest to compare the results of this three-dimensional model with the exact one-dimensional Zimm-Bragg 
model discussed in Introduction. This is done in Fig. (fTTj) for a = 10~ 3 . Note, that in the Zimm-Bragg language the e 
used here corresponds to — A/ = Ins. In the Zimm-Bragg model the increase of Or is steeper. This can be explained 
by the fact that in the three-dimensional model the entropy of the flexible and the stiff segments is explicitly taken 
into account, which hinders the generation of stiff segments. On the other hand, the mean-field character of the 
three-dimensional model might also weaken the cooperativity effect. 



C. Transition from amorphous to liquid-crystalline globule 



The main focus of this work lies on the crossover from a disordered amorphous globule with a low or moderate 
fraction of stiff segments to an ordered liquid-crystalline globule with a very high fraction of stiff segments . In Fig. (jT2J) , 
on the left, the fraction of stiff segments Or is plotted as a function of the selective two-body interaction parameter x, 
which controls the strength of the additional interaction between stiff segments and therefore models selective solvent 
conditions (hydrophobicity). On the right the nematic order parameter S is plotted as a function of \. The nematic 
order parameter is given by 



S = J d 3 r J d 2 uP 2 (cos0)^ 2 (r,u) 
d 3 r ip 2 (r) <<p (r) + V5ip 2 (r) 



(3.3) 
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FIG. 11: Or is plotted as a function of e for a — 10 3 . The continuous curve shows the numerical result for v = —0.2 and 
N = 10 . The dashed curve shows the result of the one-dimensional Zimm-Bragg theory. 



Throughout the entire section, the non-selective two-body interaction parameter is set to v — —0.2. The two plots 
in Fig. (fT2"]) demonstrate that the onset of nematic order and the increase in the fraction of stiff segments occur 
simultaneously. 
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FIG. 12: The fraction of stiff segments Or (left) and the nematic order parameter S (right) are plotted as functions of \ f° r 
a = 10~ 4 , N = 9.5 • 10 3 and e = g = 0. 



This transition occurs without an explicit angle-dependent alignment interaction, that is g — 0. The transition is 
triggered by a subtle interplay of the entropy contribution (surface energy), represented by the derivative terms in 
Eq. (|2.20p . and bulk interaction energy, represented by the x-term. This surface energy has an entropic nature since 
the conformational set of surface segments is constrained. In a simple homopolymer globule it is isotropic^. For the 
rod-coil copolymer the surface energy is anisotropic and the surface tension in p-direction is smaller than the one in 
z-direction. That is why the system tries to maximise its lateral surface in p-direction and minimize it in z-dircction, 
i.e. a nematic, cigar shaped, liquid-crystalline globule is formed. 

To demonstrate how the shape of the globule changes during the transition, Fig. (Tl3|) shows a colour-coded plot of 
the local density in g-z space. The center of the globule is in the bottom left corner of each picture, g is increasing 
from left to right and z is increasing from bottom to top. Red indicates high density and dark blue zero density. In 
Fig. (fT3| a different set of parameters is chosen: a — 10~ 4 , N = 10 5 and e = 0.1. Below, it will be discussed in detail 
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how the nature of the crossover changes with a, e and N. The values of \ m Fig. ([T3|) are chosen such that the top 
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FIG. 13: The density of the flexible segments is shown on the left and the density of the stiff segments on the right, a — 10~ 4 , 
e — 0.1 and iV = 10° for all plots, x ~ —0.0138 in the top line, % = —0.0812 in the middle line and x = —0.18 in the bottom 
line. 

two pictures (x = —0.0138) show the system at the transition point between coil and globule, the middle two pictures 
(x = —0.0812) show the system at the transition point between amorphous and liquid-crystalline globule (as defined 
below) and the bottom two pictures (x = —0.18) show the system deep in the liquid-crystalline globule phase. The 
transition from amorphous to nematic liquid-crystalline globule is a crossover of finite width and it is a rather obvious 
choice, to define the transition point as the inflection point of the S'(x)-curve in Fig. (fT2|) . At the transition point 
between coil and globule the system is spherical and has a very broad surface layer of decaying density. Although the 
density of the helical segments shown on the right is higher than the density of the flexible segments shown on the 
left, their distribution and the shape of the profile is very similar. 

At the transition point between amorphous and liquid-crystalline globule the system adopts a slightly cylindrical 
shape indicating the onset of nematic order. It can also be seen that the density maximum of the flexible segments 
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is not in the center of the globule denoting a repulsion of flexible segments from the center to the surface layer. 
The surface layer is now much narrower. Deep in the liquid-crystalline phase the globule has developed a strongly 
asymmetric cylindrical shape indicating strong nematic order. The repulsion of flexible segments from the center 
towards the surface layer can be seen clearly and the surface layer is now very narrow. 



N -dependence 



In this subsection we study how the total chain length N influences the transition from amorphous to liquid- 
crystalline globule. Fig. (fl"4f shows the fraction of stiff segments as a function of x f° r f° ur different chain lengths. 
The crossover from an amorphous globule with moderate number of stiff segments to a liquid-crystalline globule with 
very high number of a stiff segments becomes sharper with decreasing chain length, which, at first sight, is a rather 
unusual and surprising behavior. 
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FIG. 14: The fraction of stiff segments Or is plotted as a function of \ for different chain lengths, 
sharper with decreasing chain length. 



The crossover becomes 



Fig.(|15p shows the nematic order parameter 5 as a function of x f° r the f° ur chain lengths. It can be seen that the 
increase in the order parameter is stronger for shorter chains. For N = 5 • 10 5 and \ — —0.214 the system is already 
so deep in the globular phase that the corresponding value of \i, which keeps the chain length fixed at N = 5 • 10 5 , is 
very close to the value at which 2V(/z) diverges (see Fig.Q). The tuning of the chemical potential fi to ensure fixed 
chain length N becomes therefore numerically impossible for higher values of \x\- Although it might be difficult to 
see, the corresponding 6R,(x)-curve in Fig. (fl"4"|) also ends at x = —0.214. 

To demonstrate how the shape of the liquid-crystalline globule changes with chain length N, Fig. (fTo| shows colour- 
coded density plots of pc and pn for the four different values of the chain length at x = —0.18. The plots illustrate 
that the system changes from a cigar-like shape for N = 10 4 towards an almost cylindrical shape for N = 5 • 10 5 . 
Fig. (fl"6|) also demonstrates that the surface layer becomes smaller for larger systems, which indicates that a larger 
system is deeper in the globular state at the same value of x an d v but at the same time showing less nematic order. 

With increasing system size N the liquid-crystalline globule phase eventually dies out (for g = 0) because in the 
limit N — > oo the contribution of the surface terms in Eq. (|2.20p vanishes. Since for g = these terms are the only 
non-spherical symmetric terms, the system cannot adopt a liquid-crystalline globular state any more. Without explicit 
alignment interaction a liquid-crystalline globule with nematic order can therefore only form in finite systems. The 
crossover becomes sharper for smaller systems because the entropic surface terms become more important compared 
to the isotropic interaction terms. The total value of the isotropic bulk interaction energy roughly scales as the volume 
of the globule, whilst the surface energy scales as the surface area of the globule. This shows that the transition to a 
liquid-crystalline polymer globule is actually driven by entropy, due to the entropic origin of the surface energy. 
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FIG. 15: The nematic order parameter S is plotted as a function of \ f° r different chain lengths. 



Finally, in Fig. (|17|) . the average length Lr of the stiff parts is plotted as a function of x f° r N — 10 and N = 5- 10 5 . 
The figure shows that the average length of the rods is equal to 1 for very small values of x an d e = 0. For N = 10 4 
it stays almost equal to 1 until the transition point is reached and then shows a rather strong increase. For very high 
values of x it should reach an asymptotic value. However, for x > 0.25 the fraction of flexible segments (1 — Or, see 
Fig. (fT2")) ) becomes unphysically small indicating that the self-consistent field treatment reaches the limit of its validity. 
For N = 5 • 10 5 the crossover to a stronger increase in average rod length with — x is much smoother analogous to Or 
and S (see Figs. ([Til [TBI). 

The stronger increase in average rod length for shorter polymers corresponding to a sharper isotropic-nematic 
transition indicates that this transition is also a cooperative process. To study the influence of cooperativity, the 
dependence of the isotropic-nematic transition on a is discussed in the next subsection. 

2. a- dependence 

The crossover from amorphous globule to liquid-crystalline globule is a cooperative process. Fig. fTg)) shows that 
the crossover becomes sharper with increasing cooperativity (decreasing a) as expected. For very high cooperativity 
(a = 10~ 5 ) and no selective interaction (x ~ 0), the fraction of stiff segments (Or = 0.2) is much smaller than 
1/2. When the transition point is approached, the fraction of stiff segments increases rapidly. For lower cooperativity 
(larger a), not only is the crossover much smoother but the fraction of stiff segments at x ~ is much higher (Or « 0.5 
for a = 10~ 3 and a = 5 • 10~ 4 ). Both features clearly show the cooperativity of the transition. 

The behavior of the nematic order parameter S for different values of a is shown in Fig. tfT!)]) . For higher a the 
increase in nematic order is less steep and the transition point is significantly shifted to higher values of |%|. For 
a = 5 • 10~ 4 and higher, the system only develops a slight onset of order and S remains very small even deep in the 
globular state at \ — —0.25. This again demonstrates the role of cooperativity in the formation of a liquid-crystalline 
globule. 

It is therefore of interest to investigate the average rod length Lr, which should also strongly depend on the 
cooperativity. Fig. (|2"D)) shows that the increase of Lr with \x\ after the transition point becomes steeper for smaller 
values of a as one would expect. The isotropic-nematic transition is enhanced by cooperativity and even becomes 
impossible if the cooperativity is too small. As can be seen from Fig. (f2"0"|) , only a cooperative system forms long rods 
and not just many very short ones. Remember, that a is associated with the energy penalty for a boundary between 
rod and coil, as discussed in the introduction . The smaller a is, the larger is this energy penalty and the more 
favorable is the formation of long rods rather than short ones. It is intuitively clear that long rods align much easier 



17 



1 1 




latmv 

U.LII 



m-oc 

U.UJ 



ItXGGbl 

luoi6,' 

Ia00£3i 



a ct-1 1 




r;.2r;9 



FIG. 16: The density of the flexible segments is shown on the left and the density of the stiff segments on the right, a = 10 
e = = 0.0 and x = -0.18 for all plots. From top to bottom TV increases as 10 4 , 3 • 10 4 , 10 5 and 5 • 10 5 . 



than very short ones. When the average rod length stays roughly equal to 1 for high values of \x\, alignment cannot 
happen. Hence cooperativity is important to drive this transition. 
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FIG. 17: The average length Lr of the stiff parts is plotted as a function of x f° r N = 10 4 and N — 5 ■ 10 5 . 
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FIG. 18: The fraction of stiff segments is plotted as a function of \ for different values of a and N = 10 5 , e = g = 0. The 
crossover becomes sharper with decreasing a, i.e. increasing cooperativity. 



3. e-dependence 

In this subsection the e-dependence of the transition from amorphous to liquid-crystalline globule is investigated 24 . 
In Section IIB, it was shown that the fraction of stiff segments increases with increasing energy gain e per stiff segment. 
A higher value of e yields a higher offset of 0r at x = which should in turn lead to a smaller value of |x| at the 
transition point. 

Figs. pTl |22|) demonstrate this behavior. For all plots in this subsection the cooperativity parameter is set to 
fj = 10~ 4 and the total chain length to N — 10 5 . The explicit rod-rod alignment interaction is still switched off 
(g = 0). The onset of the transition is indeed shifted to lower values of |x| with increasing energy gain per helical 
segments e. This is due to an increase of bulk interaction energy for fixed x with increasing number of stiff segments. 
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FIG. 19: The nematic order parameter S is plotted as a function of x f° r different values of a and N = 10 s . The transition 
point is shifted to higher values of \x\ with increasing a. 




The curves for e > 0.1 in Figs. (|2"Tl |2"2")) start at non-zero values of \. These values of \ correspond to the transition 
point between open chain regime and globular regime, similar to the discussion in Section IIA. The main difference 
in this case is that at v = —0.2 the non-selective attractive two-body interaction is not strong enough to drive the 
system into the globular state. The additional selective two-body interaction between the stiff segments only (with 
interaction parameter \) is needed to drive the system into the globular state. That this is only the case for e > 0.1 
can be explained as follows (see also Section IIB). If there is no additional selective interaction energy which favors a 
compactification of the stiff segments, the stiffening of parts of the chain due to an increase of 0r with increasing e 
pushes the chain segments further apart from each other and therefore leads to a more open structure. For e > 0.1, 
the system is thus pushed into the open chain regime at x — 0- On the other hand, for e > 0.17 the system can be 
driven directly from open chain to liquid-crystalline globule as soon as \x\ exceeds some crossover value. 

These findings permit the computation of a complete phase diagram of the rod-coil copolymer in e~x space, see 
Fig. ([23)) . But before this is done, it has to be checked whether the definition of the transition point between coil 
and globule as the minimum of the N (/i)-curve is also reasonable in the case of strong nematic order. In Section 
IIA the radius R of the globule (defined as the point at which the total density has decreased to p{R) = 10~ 3 po) 
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FIG. 21: Fraction of helical segments 0r as a function of \ for different values of the energy gain per helical segment e. 
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FIG. 22: Nematic order parameter S as a function of x f° r different values of e. 



was plotted as a function of v, see Fig. ([7]). It showed a rapid increase when the transition point between globule 
and coil was approached. Here the globule shows nematic order and has an asymmetric shape, i.e. it is necessary 
to distinguish between g- and z-direction. The extensions of the globule in p-direction R g and in z-direction R z are 
defined - analogously to R in the case of a spherical globule - as p(R e ,0) — p(0,R z ) — 10~ 3 po. In Fig. (f23|) R e and 
R z are plotted as functions of x f° r e — 0.25. For this choice of e, the system already shows strong nematic order 
at x = —0.0396 which corresponds to the transition point between globule and coil. Fig. (|2"3"|) demonstrates that the 
definition of the transition point between open chain and globular regime as the minimum of the N-p curve is valid 
even if the system shows strong nematic order. Both curves in Fig. (|23[) show a rapid increase when approaching the 
transition point from the right. The different values of R e and R z reflect the cylindrical shape of the globule due to 
the nematic order of the stiff segments. For e = 0.25 no isotropic globular phase exists and the crossover leads directly 
to a liquid-crystalline globule. 
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FIG. 23: The extensions of the globule in p-direction R e and in z-direction R z are plotted as functions of x for e = 0.25. 



For e — 0.1 on the other hand, the system shows no nematic order at the transition point between open chain and 
globule (x — —0.0138). The system therefore undergoes two transitions. First from an open chain to an amorphous 
globule and then at higher \\\ from an amorphous to a liquid-crystalline globule. In Fig. (|24[) R g and R z are plotted 
as functions of x- When the transition point between open chain and globule (at x — —0.0138) is approached from 
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FIG. 24: The extensions of the globule in ^-direction R e and in z-direction R z are plotted as functions of \ for e = 0.1. 

the right, R B and R z show the expected strong increase. When the transition point between amorphous and liquid- 
crystalline globule (at x — —0.0812) is approached from the left, R B and R z show a different behavior. In a small 
interval (corresponding to the width of the transition) R z increases whilst R g decreases even stronger than before. In 
this interval the asymmetric shape of the globule is developed. This behavior will be discussed further in Subsection 
D. 

As already mentioned above, it is now possible to compute a complete phase diagram in e-x space. This phase 
diagram is shown in Fig. (|25[) . The triangles are the transition points between open chain and globule. The squares 
are the transition points between amorphous and liquid-crystalline globule. Note, that the points plotted in the phase 
diagram, Fig. (|2"5|) , are what is defined above as the points of rather broad crossover regions. Therefore the boundaries 
in the phase diagram have to be understood as center lines of broader regions in which the crossover from one phase 
to the other occurs. The little arrows in Fig. ((25|) correspond each to one of the rows of the color coded density plots 
shown in Fig. (fT3"|) , helping to illustrate how the density profiles look at these points in the phase diagram. 

4- g-dependence 

The transition from amorphous to liquid-crystalline globule occurs without an explicit angle-dependent alignment 
interaction between the rods. It is nevertheless interesting to switch the explicit rod-rod alignment interaction on 
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FIG. 25: Phase diagram of a rod-coil copolymer in e-\ space. The upper left area corresponds to an open chain, the lower left 
area to an amorphous globule and the right area to a nematic liquid-crystalline-globule. The little arrow to the left indicates 
the point in the phase diagram which corresponds to the top two pictures in Fig. {13}. The arrow in the middle corresponds to 



the middle two pictures and the arrow to the right to the bottom ones. Parameters: N = 10 , a = 10 



0. 



and study its influence on the transition. The alignment interaction is of the Maier-Saupe form and its strength is 
controlled by the interaction parameter g, see Eq. (|2.20p . It is chosen to be attractive to favor alignment of the rods. 

In Fig. (|2l)|) the fraction of stiff segments is plotted as a function of \ f° r different values of g and a = 10~ 4 , e = 0, 
N = 9.5 • 10 3 . The corresponding nematic order parameter S is shown in Fig. (|27[) . The transition is shifted to 
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FIG. 26: Fraction of helical segments Or as a function of \ f° r different values of g. 
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lower values of values of \\\ with increasing \g\. This is not suprising, since for g < the alignment term provides 
an additional incentive for the system to generate nematic order and stiff segments. But the fact that almost up 
to the respective transition points the 0ft(x)-curves lie perfectly on top of each other is quite remarkable. This can 
be explained as follows. For tp2 = the alignment term in Eq. (|2.20p is equal to zero. Therefore the system has to 
generate a finite tp2 (i.e. at least small nematic order) before the alignment interaction can have an effect on the 
system. As long as there is no nematic order in the system (ip2 = 0) the alignment term is zero and the system 
behaves as the one with g — 0. As expected, the attractive alignment interaction enhances the generation of nematic 
order and thus shifts the transition point to lower values of |x|- 



D. Free energy 

In this subsection the behavior of the effective free energy in the crossover region is investigated. The effective 
saddle point grand potential is given by Eq. (|2.20p . With the definition of the chemical potential used here the 
corresponding effective saddle point free energy is given by 

F = fl-fxN. (3.4) 

In Fig. (|28|) the free energy F is plotted as a function of \ for a = 10~ 4 , TV = 9.5 • 10 3 and e = g = 0. These parameters 
are the same as the ones used in Fig. (TT2"|) . Fig. (|2"8|) demonstrates that the free energy is changing its slope in the 
crossover region around the transition point as would be expected for a crossover. At the beginning of Section IIC 
the occurrence of the transition from amorphous to liquid-crystalline globule was explained in terms of the interplay 
between surface energy and bulk interaction energy. Because of the anisotropy of the surface contributions which 
originate from the entropy of the rods, the globule tries to minimise its surface in z-direction and to maximise it in 
p-direction. Before the transition the amorphous globule has a spherical shape. The surface energies in x-, y- and 
z-direction should therefore be all the same. For the cylindrical coordinates used here that implies F^ ur! = 2F™ l{ . 
Fig. ([29f shows that this is indeed the case. But also in the liquid-crystalline globule regime after the transition their 
ratio is roughly equal to 2 as can be seen from Fig. (j2"9"|) . The extension of the globule in g- and z-direction on the 
other hand is now very different. It is therefore instructive to plot F™ and F™ normalized by the corresponding 
cross sections of the globule, which should give a feeling how the surface tension behaves. As an approximation of 
the cross section in z-direction Ft? is chosen. The cross section in p-direction is approximated by R e R z , where R e 
and R z are defined as in Subsection IIC3. In Fie. (|30|) the normalized surface contributions are plotted on the left 
hand side and R e and R z one the right hand side. The plot on the left shows that the surface energy per area in 
z-direction becomes indeed larger than the one in ^-direction in the interval in which the crossover from amorphous 
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globule to liquid-crystalline globule occurs. The corresponding extensions of the globule R e and R z plotted on the 
right illustrate the enlargement of the globule in z-dircction and the diminution in p-direction in the crossover interval. 
When the final shape of the liquid-crystalline globule has developed the curves become parallel again and decrease 
both indicating the further compactification of the entire globule. The investigations of the surface contributions to 
the free energy further visualize that it is the anisotropy of the entropic surface energy that drives the system into 
the ncmatic state. 



IV. CONCLUSIONS 



The numerical solutions of the self-consistent field equations show that the rod-coil copolymer with variable com- 
position can form three phase states, open chain, amorphous globule and nematic liquid-crystalline globule with high 
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FIG. 30: F s g url /(R e R z ) and F! nrC /(R 2 ) plotted as functions of X for a = 1(T 4 , N = 9.5 ■ 10 3 and e = g = 0. 



fraction of stiff segments. The transition between the first two states is similar to the coil-globule transition of a 
homopolymer. The formation of a liquid-crystalline globular state without explicit alignment interaction between the 
rods is a novel result and deserves further discussion. 

The formation of a liquid-crystalline globule from a rod-coil multiblock copolymer with hxed composition has been 
discussed for the first time in an early theoretical work^l. Phase states similar to the ones summarized in the phase 
diagram in Fig. (|25|) have been seen in^, where the authors consider a homopolymer in which each monomer carries 
a dipole moment and take into account explicit dipole-dipole interactions. However, the schematic phase diagram 
ir*2£ has been considered within the so-called volume approximation where the contribution of the surface energy can 
be neglected (at N —* oo)22 and the transition into the anisotropic globular state with nematic order is driven by 
dipole-dipole interactions. 

Recently, in the paper by Marenduzzo et al^ a new class of models for chain molecules has been considered . This 
model can be viewed as an elastic tube and the concept of the chain thickness (" thick polymer" ) has been introduced 
by means of a specially prescribed three-body interaction potential. Nevertheless, we should stress that on the mean- 
field level the consideration (within the so-called " chain of coins" model) is still limited to the volume approximation, 
and the explicit alignment interaction in form of an Onsager second virial term is present (see the Eq.(14) in ref. — ). 
In contrast to this in our case the formation of a LC-globule occurs without explicit alignment interactions between 
the helical parts. It is the entropic surface tension anisotropy which drives the globule in the nematic LC-state. 

In an early work^P. a cooperative helix-coil liquid-crystal transition has been found, very similar to the transition 
from amorphous to liquid-crystalline globule discussed here. The formation of nematic order is also accompanied by 
a strong increase in fraction of stiff (or helical) segments. The main difference is again that the transition is driven 
by an explicit alignment interaction of Maier-Saupe type similar to the additional interaction (g < 0) considered in 
Section IIC4 and not by the entropic surface tension anisotropy. 

Monte Carlo simulations of a simple homopolymer model on the lattice have been carried out to study the interplay 
between coil-globule transition and a-helix formation^. The phase diagram has been presented in terms of two 
energetic parameters that characterize the hydrophobic attraction between monomers that are far apart in sequence 
and the local helix-stabilizing interaction. This is reminiscent of our x — e diagram given in Fig|25l The number of 
helical units, (2), (helicity) and the number of pairwise contacts, (NC), are used as criteria to distinguish among 
different phases. Following these criteria the authors are able to see four different phases. Nevertheless it is difficult 
to compare these results with our findings because the order parameters we use are substantially different form the 
ones used in reff« . 

In the context of secondary and tertiary structure formation in proteins the interplay of helix formation and liquid- 
crystalline order has been studied, see for instanc o 31 i 32 . It was shown that liquid-crystalline ordering enhances the 
number of helical segments as well as the average length of the helices which is qualitatively similar to the results of 
the model presented here. This indicates that the model of a rod-coil multiblock copolymer with variable composition 
might be a good candidate to give a simple explanation for the formation of helix bundles in certain globular proteins. 
Both, simulations^ 3 - and experiments 34 show that proteins can adopt not only the native and completely denatured 
state (open chain) but also so-called premolten and molten globular states. For helix-bundle proteins the premolten 
globule, which does not show any order of the helices, corresponds to the amorphous globule in our model. The molten 
globule with ordered helices but without native contacts (and therefore also without the characteristic helix-helix angle 
found in the native state) corresponds to the liquid-crystalline globule. Our model also shows that during the transition 
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to a liquid-crystalline globule not only the amount of helical segments increases strongly but the globule also becomes 
more compact. The experimentally observed^ correlation between the amount of secondary structure elements and 
compactness of proteins mentioned in the introduction might therefore also be explained by this transition, at least 
in the case of helix bundle proteins. 

Irrespective of the possible application to helix bundle formation in proteins, the model which we have discussed in 
this paper provides a relatively simple example of the general interplay between secondary structure (helices or stiff 
rods) and tertiary structure (liquid-crystalline order) in homopolypeptides. 
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